import matplotlib.pyplot as plt
import numpy as np
from scipy.stats import binom

# 设定二项分布的两个参数
n = 100       # 实验次数
p = 0.02      # 成功概率

# 计算二项分布的概率密度函数
x = np.arange(0, 11)             #做等差数列，直观上是坐标x轴
y = binom.pmf(x, n, p)           #二项分布的概率密度函数，内置公式为书上分布律。

# 绘制二项分布图
plt.scatter(x, y)                #散点图
plt.plot(x, y, linestyle='solid',color='blue', label='Probability Density')    #点之间用线段连接，增加图例
plt.xlabel('Successes')
plt.ylabel('Probability Density')
plt.title('Binomial Distribution')
plt.show()